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We analyze the evolution equation describing all multiple hard pomeron exchanges in 
a hadronic or nuclear structure functions that was proposed earlier in [Q]. We construct 
a perturbation series providing us with an exact solution to the equation outside of the 
saturation region. The series demonstrates how at moderately high energies the corrections 
to the single BFKL pomeron exchange contribution which are due to the multiple pomeron 
exchanges start unitarizing total deep inelastic scattering cross section. We show that as 
energy increases the scattering cross section of the quark-antiquark pair of a fixed transverse 
separation on a hadron or nucleus given by the solution of our equation inside of the satura- 
tion region unitarizes and becomes independent of energy. The corresponding F2 structure 
function also unitarizes and becomes linearly proportional to Ins. We also discuss possible 
applications of the developed technique to diffraction. 



I. INTRODUCTION 



The Balitsky, Fadin, Kuraev and Lipatov (BFKL) equation resums all leading logarithms of Bjorken 
X for hadronic cross sections and structure functions. The solution of the BFKL equation grows like a power 
of center of mass energy s, therefore violating the unitarity bound at very high energies [^,|[. This is one 
of the major problems of small x physics, since the Froissart bound [^jj^ states that the total cross section 
should not raise faster than In^ s at asymptotically high energies. There is a general belief that the unitarity 
problem could be cured by resumming all multiple BFKL pomeron exchanges in the total cross section, or, 
equivalently, in the structure function. 

As was argued by Mueller in |6| multiple pomeron exchanges become important at the values of rapidity 
of the order of 
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with a the strong couphng constant, which is assumed to be small, and ap — 1 



In 2. This result could 



be obtained if one notes that one pomcron contribution to the total cross section of, for instance, onium- 
onium scattering, is parametrically of the order of a^exp[(ap — 1)^^] and the contribution of the double 
pomeron exchange is Q;^exp[2(Q;p — 1)1"]. Since multiple pomeron exchanges become important when the 
single and double pomeron exchange contributions become comparable, we recover Eq. (|l|) by just equating 
the two expressions. 

After completion of the calculation of the next -to-leading order corrections to the kernel of the BFKL 
equation (NLO BFKL) by Fadin and Lipatov Q, and, independently, by Camici and Ciafaloni it was 
shown in B] that due to the running coupling effects these corrections become important at the rapidities 
of the order of 
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One can see that Yu <ti Ynlq foi' parametrically small a. That implies that the center of mass energy at which 
the multiple pomeron exchanges become important is much smaller, and, therefore, is easier to achieve, than 
the energy at which NLO corrections start playing an important role. That is multiple pomeron exchanges 
are probably more relevant than NLO BFKL to the description of current experiments jio). Also multiple 
pomeron exchanges are much more likely to unitarize the total hadronic cross section. In this paper we are 
going to propose a solution to the problem of resummation of multiple pomeron exchanges and show how 
the BFKL pomeron unitarizes. 




FIG. 1. An example of the pomeron "fan" diagram. 

In the previous paper on the subject we proposed an equation which resums all multiple BFKL pomeron 
1^,^ exchanges in the nuclear structure function. In deriving the equation we used the techniques of Mueller's 
dipole model As was proven some time ago the dipole model can provide us with the 

equation which is exactly equivalent to the BFKL equation In [|l| it was shown that the dipole 

model methods, when applied to hadronic or nuclear structure functions, yield us with the equation which 
in the traditional Feynman diagram language resums the so-called pomeron "fan" diagrams in the leading 
logarithmic approximation, an example of which is depicted in Fig. y for the case of deep inelastic scattering. 

To review the results of ^ let us consider a deep inelastic scattering (DIS) of a virtual photon on a hadron 
or a nucleus. An incoming photon splits into a quark-antiquark pair and then the qq pair rescatters on the 
target hadron. In the rest frame of the hadron all QCD evolution should be included in the wave function 
of the incoming photon. That way the incoming photon develops a cascade of gluons, which then scatter 
on the hadron at rest. In [|l| this gluon cascade was taken in the leading longitudinal logarithmic (Inl/x) 
approximation and in the large Nc limit. This is exactly the type of cascade described by Mueller's dipole 
model pl]-p^ . In the large Nc limit the incoming gluon develops a system of color dipoles and each of them 
independently rescatters on the hadron (nucleus) M . The forward amplitude of the process is shown in Fig. 
^. The double lines in Fig. |^ correspond to gluons in large Nc approximation being represented as consisting 
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of a quark and an antiquark of different colors. The color dipoles are formed by a quark from one gluon and 
an antiquark from another gluon. In Fig. || each dipole, that was developed through the QCD evolution, 
later interacts with the nucleus by a series of Glauber-type multiple rescatterings on the nucleons. The 
assumption about the type of interaction of the dipoles with the hadron or nucleus is not important for the 
evolution. It could also be just two gluon exchanges. The important assumption is that each dipole interacts 
with the target independently of the other dipoles, which is done in the spirit of the large Nc limit. 

For the case of a large nucleus independent dipole interactions with the target nucleus are enhanced by 
some factors of the atomic number A of the nucleus compared to the case when several dipoles interact with 
the same nucleon This allowed us to assume that the dipoles interact with the nucleus independently. 
The summation of these contribution corresponds to summation of the pomeron "fan" diagrams of Fig. 
|l|. However, in general there is another class of diagrams, which we will call pomeron loop diagrams. In a 
pomeron loop diagram a pomeron first splits into two pomerons, just like in the fan diagram, but then the two 
pomerons merge back into one pomeron. In other words the graphs containing pomerons not only splitting, 
but also merging together will be referred to as pomeron loop diagrams. In a large nucleus these pomeron 
loop diagrams could be considered small, since they would be suppressed by powers of A. This approximation 
strictly speaking is not valid for the calculation of a hadronic wave function. At very high energies it also 
becomes not very well justified even in the nuclear case for the following reason. The contribution of each 
additional pomeron in a fan diagram on a nucleus is parametrically of the order of a^A^/'^e^"^"^^^, which 
is large at the rapidities of the order of Yu given by Eq. (|^). The contribution of an extra pomeron in a 
pomeron loop diagram is of the order of a^e'^"^"^''^, which is not enhanced by powers of A and is, therefore, 
suppressed in the nuclear case. However, one can easily see that at extremely high energies, when a^e*^"^"^'^ 
becomes greater than or of the order of one, pomeron loop diagrams become important even for a nucleus, 
still being smaller than the fan diagrams. 




FIG. 2. Dipole evolution in the deep inelastic scattering process as pictured in Ref. [1]. The incoming virtual photon 
develops a system of color dipoles, each of which rescatters on the nucleus (only two are shown). 

Resummation of pomeron loop diagrams in a dipole wave function seems to be a very hard technical 
problem, possibly involving NLO BFKL or, even, next-to-next-to-leading order BFKL kernel calculations 
. Nevertheless, it is the author's belief that these pomeron loop diagrams are not going to significantly 
change the high energy behavior of the structure functions. Still in the rigorous sense the results of this 
paper should be considered as resumming the fan diagrams only, which is strictly justified only for a large 
nucleus. 

In 0, rewriting the F2 structure function of the target as a convolution of the wave function of the qq 
fiuctuations of a virtual photon with the propagator of the quark-antiquark pair through the nucleus, we 
obtained 

F2(x,Q2) = / ^!^ci,(xoi,z) d% 7V(xoi,bo,y), (3) 

An^aEM J 27r 

where the photon's wave function 

$(xoi,z) = $t(xoi,z) + $l(xoi,0) (4) 

consists of transverse 
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and longitudinal 



$z.(xoi, z) = ^^^^ AQ'z\l zf KUxoia) 



(5b) 



contributions, with = Q^z{l ~ z), and we assumed that quarks are massless and have only one flavor 
and unit electric charge. As always is the virtuality of the photon and x is the Bjorken a;-variable. z 
is the fraction of the photons light cone momentum carried by the quark, Xqi = Xi — Xq is the transverse 
separation of the quark-antiquark pair, bo is the impact parameter of the original pair and Y ~ ln(l/x) 
is the rapidity variable. iV(xoi,bo,F) is the propagator of the quark-antiquark pair through the nucleus, 
or, more correctly, the forward scattering amplitude of the qq pair with the nucleus (hadron). It obeys the 
equation j|] 



iV(xoi,bo,F) = -7(xoi,bo) exp 



AaCp , / xoi 
m 



Y 
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dy exp 



f x'^ 1 1 1 

X / d'^x2 [2jV(xo2,bo + -xi2,y)-A^(xo2,bo + -xi2,2/)iV(xi2,bo + -Xo2,y)], (6) 

with 7(xoi, bo) the propagator of a dipole of size Xoi at the impact parameter bo through the target nucleus 
or hadron, which was taken to be of Glauber form in : 



7(xoi,bo) = exp 



-xl^AxG{x,l/x.l^) 



1. (7) 



here S± is the transverse area of the hadron or nucleus, A is the atomic number of the nucleus, and xG is 
the gluon distribution of the nucleons in the nucleus, which was taken at the two gluon (lowest in a) level 
[p7| . Eq. (^ is written for a cylindrical hadron or nucleus. In the large Nc limit Cp ~ Nc/2. In Eq. (||) p 
is an ultraviolet regulator 

Eq. (||) resums all multiple hard ponieron exchanges in the F2 structure function. The linear in N term 
on the right hand side of Eq. (^ corresponds to the usual BFKL evolution , whereas the quadratic term 
in A'^ is in a certain way responsible for triple pomeron vertices. This quadratic term on the right hand side 
of Eq. (^ introduces some damping effects, which at the end unitarize the growth of the BFKL pomeron. 
We note that an equation similar to (^ was suggested by Gribov, Levin and Ryskin (GLR) in momentum 
space representation in jl^ and was subsequently studied by Bartels and Levin in However, usually 

by the GLR equation one implies an equation summing all fan diagrams in the double leading logarithmic 
approximation. Eq. (^) resums all the fan diagram in the leading longitudinal logarithmic approximation. In 
the double logarithmic limit (large Q^) Eq. (^ reduces to the GLR equation. We also note that an equation 
very similar to Eq. ^ has been proposed by Balitsky in to be used in describing the wavefunctions of 
each of the quarkonia in onium-onium scattering, similarly to |ll2| , pl| . A renormalization group techniques 
based on McLerran-Venugopalan model have been developed by Jalilian-Marian, Kovner, Leonidov, 
and Weigert leading to an equation which is also supposed to resum multiple hard pomeron exchanges 
psf . Unfortunately the resulting equation |Q differs from our Eq. (||) . An equation summing all multiple 
pomeron diagrams in the double logarithmic limit has been proposed by Ayala, Gay Ducati and Levin |24[| , 
which is also different from the appropriate limit of our Eq. (^) . 

In this paper we are going to analyze Eq. (^. In Sect. II we will find a solution of Eq. (||) in the form 
of a series in Eq. ( ^ ) which will be constructed by solving Eq. (^) perturbatively in the momentum space. 
The details of transforming Eq. into momentum space will be given in the Appendix A with the result 
given by Eq. The terms of the series will correspond to multiple pomeron exchange contributions. 

Partial sums of the series will be plotted in Fig. |3[ which shows how multiple pomeron exchanges unitarize 
the amplitude. We will argue that the kinematic region for which the series is convergent corresponds to the 
case when the DIS cross sections and structure functions have not reached saturation. In a certain sense we 
will define the saturation region as the kinematic regime where the series ceases to converge. By resumming 
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the terms of the series corresponding to exchanges of large numbers of pomerons we construct a function 
which we can analytically continue outside the region of convergence of the series of Eq. (|3^) . That would 
provide us with the ansatz for the solution of Eq. ( |l0|) outside of the saturation region [see Eq. (^)] . We 
will prove that this ansatz actually is the only possible energy independent solution of Eq. (p^. 

To transform our results into coordinate space on would have to separate the momentum integration 
in two regions, one of them being above and the other one being below the saturation scale. Since we 
will find only the asymptotic solution inside the saturation region our results may not be very precise in 
determining the coordinate space amplitude. A numerical solution of Eq. would probably allow one to 
achieve better precision. Nevertheless for asymptotically high energies we can transform the amplitude into 
coordinate space obtaining Eq. (52). Thus we will show that multiple pomeron exchanges unitarize the 
forward amplitude, making it independent of energy at high energies, which is shown in Fig. ^. 

In Sect. Ill we will analyze the high energy asymptotics of the F2 structure function. We will show 
that, contrary to some expectations, F2 structure function does not become independent of e nerg y at high 
energies. It becomes linearly proportional to Ins, still satisfying the Froissart bound [see Eq. ip(\)]. 

In Sect. IV we will discuss how in the framework of the dipole model one can describe diffractive scattering 
by including the multiple pomeron evolution in rapidity of Eq. (^) in the quasi-elastic structure functions 
considered in |p5|. W e will derive an expression for the quasi-elastic structure function (|4|) and discuss 
its asymptotics(|66|) . Finally, we will conclude in Sect. V by summarizing our results and discussing the 
possible implications of the developed techniques on different issues of small- a; physics. 



II. PERTURBATIVE SOLUTION OF THE EVOLUTION EQUATION 

We start by considering a DIS on a large cylindrical nucleus. This assumption about the size and geometry 
of the nucleus is being done for simplicity of calculations and our results could be easily generalized to a large 
nucleus of any shape. If the nucleus is sufficiently large we can assume that the forward amplitude of the 
scattering of the original quark-antiquark pair on the nucleus, N{xoi, bp, Y), is a slowly varying function of 
the impact parameter of the virtual photon bp. This is equivalent to assuming that the transverse extent of 
the dipole wave function is much smaller than the size of the nucleus, and that the impact parameter between 
the virtual photon and the nucleus is smaller than the radius of the nucleus R. We note that in order to 
construct iV(xoi, bo, Y) in [|l[ we integrated over the impact parameters b^'s of the nucleons with which each 
of the pomerons in Fig. [I] interacted. That way bo does not have a meaning of the impact parameter of a 
pomeron exchange anymore and bo dependence in iV(xoi,bo,y) becomes purely "geometrical" . In a more 
formal language we assume that x± <C -R, where xj_ is the typical transverse size of the dipole wave function, 
and bo < R [bo is zero at the center of the nucleus). If these conditions are satisfied we may neglect X12 and 
X02 compared to bo on the right hand side of Eq. (^, which means that the changes in N{xqi, bo, Y) when 
bo is varied by xj^ are negligibly small. Thus suppressing the bo dependence we write 7V(xoi, Y) instead of 
N{xoi,ho,Y) in Eq.(||). Then differentiating Eq. (^) with respect to Y we can rewrite it as 



a7V(xoi,y) _ 2aCF 
dY ^ 7r2 



2^02'^12 



27r^^(xoi-Xo2)lnP^ 



iV(xo2,F) 



d2x2^^iV(xo2,r)7V(xi2,r), (8) 

^02^^12 

with iV(xoi, 0) = — 7(xoi) as the initial condition, where we also suppressed (but did not neglect) the impact 
parameter dependence of 7. Now one can see that the evolution equation (|^) itself in the large nucleus 
approximation does not impose any impact parameter bo dependence. The latter is only included in the 
initial conditions given by 7(xoi,bo). Since 7 is a slowly varying function of bo (constant for a cylindrical 
nucleus) our assumption about the slow dependence of N on bo is self consistent. Physically this implies 
that before the collision, when the incoming virtual photon develops its dipole evolved wave function, it does 
not have any information about the nucleus. The information about the thickness of the nuclear medium as 
a function of the impact parameter comes through 7 when the dipoles interact with the nucleus. 
Let us rewrite 




5 



N{x±,Y) = x'i J —e'^-''N{k,Y)^xl J dk k Jo{kxj_) N{k,Y). (9a) 
We note that the inverse of this transformation is 

N(k, Y) = / e-'"^-^ N{x^,Y) = / Mkx^) N{x^,Y). (9b) 

J Znx^ Jq x± 

Since N{x±,Y) does not depend on the direction of x^, N{k,Y) does not depend on the direction of the 
relative transverse momentum of the qq pair k^. That aUowed us to simphfy the integrations in Eqs. ( pa| ) 



and (9b). Now Eq. (g|) becomes 



dN{k,Y) _ 2aNc ( d \ aN^ 



xkinrTpv(fc,r)-— iiv^(fc,r), (lo) 



dY -n \ 91nfc 

where 

x(A)=^(i)-i^(i-^)4v.Q) (11) 

is the eigenvalue of the BFKL kernel |,|,|ll| with V(A) = r'(A)/r(A). In Eq. 1^ the function x(A) is taken 
as a differential operator with A = —9/9 In /c acting on N(k^Y\ We also put iVc/2 instead of Cp everywhere 
in the spirit of the large iV^ approximation. The details of obtaining Eq. ( p^ from Eq. (H) are presented in 
the Appendix A. Eq. ( p^ ) shows explicitly that the first, linear, term on its right hand side gives the usual 
BFKL equation, whereas the quadratic term corresponds to the triple pomeron vertex, which unitarizes the 
growth of the BFKL pomeron. 



One can easily prove that Eq. (10) has a unique solution satisfying given initial condition at Y — Q. The 
strategy of the proof is conventional: assume that there are two different solutions of Eq. ( |lo|) satisfying 
given initial conditions. Then one can easily derive an equation for the function given by the difference of 
these solutions with the initial condition stating that the function is equal to zero when y = 0. After a little 
analysis of the resulting equation with that initial condition one could easily see that the difference between 
the two solutions of the original Eq. (^0|) is zero. That way one proves uniqueness of the solution of Eq. 

We want to find a solution of Eq. (|10|) satisfying the initial condition given by the BFKL pomeron 
contribution at relatively "small" rapidity Y ~ 1/a. Unfortunately finding an exact analytical solution 
of Eq. ( p^ seems to be a very difficult task, partly because Eq. ( |lO| ) is non-linear, partly because of a 
complicated structure of the BFKL kernel x(A). Instead we are going to construct a series resulting from 



the perturbative solution of Eq. (10), which, in the region where it is convergent, provides us with the exact 
solution of Eq. (|o|). As we will show below the region of convergence of that series corresponds to the 
regime when the DIS cross sections and structure functions are not saturated. In the saturation regime the 
series diverges, but allows us to construct an asymptotic solution by analytical continuation. 

Let us start by assuming that at Y ^ 1/a the function N is very small, iV <C 1. This allows us to neglect 
the quadratic term in Eq. (pO|), since it would be of higher order in N . Then Eq. ( p^ ) becomes 

^*''^-^^"'-^""'xf-^U.(M-). (12) 



dY \ ainfc, 

which corresponds to the usual BFKL equation . The solution of Eq. (|l^) is given by 



7Vi(fc,r) = exp 



2aN^ ( d 
Yx 



dink 



C'ik), (13) 



with C{k) being some unknown function of k, to be specified by initial conditions. 
Rewrite C{k) as a Mellin transformation 



dX / fc " ^ 
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where A ^ Aqcd is some scale characterizing the nucleus or hadron. Here the integration in A runs parallel 
to the imaginary axis to the right of all the singularities of C\. Substituting Eq. (|lj) into Eq. ( p^ ) we 
obtain 



Ni(k,Y) = 



2ni 



exp 



2aNr 



Cx. 



(15) 



Assuming that the transverse momentum of the qq pair k is not very large, such that ln(fc/A) ^ oYNc 
and that Cx{Y) is a slowly varying function of A we come to the conclusion that the integral in Eq. ( p^ 
is dominated by the saddle point in the exponential in the vicinity of A = —1 in x(~A). Expanding 
;^(— A) w 21n2 + |C(3)(A + 1)^, with C(z) the Riemann zeta function, we can perform a saddle point 
approximation in the integral of Eq. (nS) . The saddle point is given by 



where we have defined 



Asp = -1 - a{k,Y), 
7rln(fc/A) 



a{k,Y) 



7aiV,C(3)y' 



which is small. The result of the integration yields 

A exp[(ap - l)Y] 



Ni{k,Y)^C, 



k ^UaNcC{3)Y 



exp 



UaNrXi3)Y 



■In' 



A 



Pi(fc,y), 



(16) 



(17) 



(18) 



a factor usually associated with the zero momentum transfer single hard pomeron exchange, which we 
denoted by Pi. C_i is determined from initial conditions. If the initial conditions are given by the two 
gluon exchange approximation, then C_i ^ a^. That way, parametrically, when Y ^ l/a, the amplitude 
iVi ~ <C 1, and our assumption of the smallness of TV at y ~ 1/a is justified. 

Eq. ( p^ ) provides us with the usual single BFKL pomeron solution. Now we are going to find corrections 
to this expression resulting from Eq. ([ic|), which will form a series 



7V(fc,y) = ^iV„(fc,F), 



(19) 



with Nn being of the order of iVf. The first correction is N2- Rewriting N = Ni + N2, substituting it back 
into Eq. (|l^), neglecting non-linear terms in N2 and employing Eq. (|lj) we obtain 



dN2ik,Y) _ 2aN^ 



dY 



d 



d\nk 



N2(k,Y) 



^N.ik,Y)^, 



(20) 



where iVi is given by Eq. ([l5|). One can easily see that the solution of Eq. ( po| ) with the initial condition 
that it disappears at small Y is 



N2ik,Y) 



dy exp 



2aNr 



d 



TT JO 

Plugging Ni of Eq. ( p^ into Eq. (|2^) and integrating over y yields 

A1+A2 



dink 



N2(k,Y) = - 



d\i d\2 ^ ^ ( k 



[Ni{k,y? 



1 



2[x(-Ai)+x(-A2)-x(-Ai-A2) 



(21) 



exp 



>^[x("Ai)+x(-A2)] 



exp( 2^yx(-Ai-A2) 



(22) 
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Eq. (20) corresponds to putting one triple pomeron vertex (given by quadratic term in TV) into the DIS 
diagram, so that N2 consists of one pomeron splitting into two. The integration over y in Eq. (|2l] ) corre- 
sponds to integration over all possible rapidity positions of the triple pomeron vertex between the inco ming 
quark-antiquark pair and the target. Thus the first exponential term in the square brackets in Eq. ( P2| ) 
corresponds to the case when the pomeron splitting occurs very close (in rapidity) to the qq pair, so that 
most of the rapidity interval is covered by two pomeron contribution. The second exponential term in Eq. 
( p^ ) corresponds to the situation when the splitting happens very close to the target, so that most of the 
interaction corresponds to the single pomeron exchange. Therefore the first exponential term in Eq. ( p^ ) 
corresponds to two pomeron contribution and the second term corresponds to the one pomeron contribu- 
tion. Similar conclusions have been reached by Navelet and Peschanski in pl[ | when analyzing two pomeron 
contribution to onium-onium scattering in dipolc model [p^ . As we will see below the two pomeron con- 
tribution gives a term proportional to exp[2(ap — 1)Y], while the one pomeron exchange is proportional to 
exp[(Q;p — 1)Y]. That way the two pomeron term dominates in A'2 and the single pomeron term could be 
safely neglected. One should also remember that A'2 is proportional to the square of C\, so that the single 
pomeron exchange term from N2 represents just an order (i.e., next-to-next-to- leading order) correction 



to iVi. Neglecting this contribution we rewrite Eq. (E2^ as 



N2{k,Y) = - -^Cx.Cx, -r 



\AJ 2[x(-Ai)+x(-A2)-x(-Ai-A2)] 



X exp 



i^[x(-Ai)+x(-A2) 



(23) 



Here again we can perform the integration over Ai and A2 in the saddle point approximation. The saddle 
points of Ai and A2 integrations would have been independent if it were not for the x(~Ai — A2) in the 
denominator. The problem is that x(A) has a pole at A = 2, and since the saddle points of Ai and A2 
integrations are located in the vicinity of —1 [see Eq. (|l^)] this could influence the exact location of the saddle 
points of the integral. However, as one can see analyzing the saddle points of Eq. (|2^), if ln^(fc/A) > aYNc 
the effect of x(— Ai — A2) could be safely neglected. This imposes an additional constraint on the possible 
values of transverse momentum k. Remember that in order to perform the saddle point integration in Eq. 
( p^ ) we had to assume that ln(fc/A) ^ aYNc- These two conditions specify the range of k in which our 
analysis is applicable. The condition ln^(A:/A) > aYNc is not very hard to satisfy, since, after all, the 
transverse momentum k should be large enough for perturbative QCD to be applicable, i.e. , ln(A:/A) ^ 1. 
We will return to this question below. 

Performing the saddle point integration in Eq. (G3) in the spirit of the above discussion one obtains 



N2{k,Y) 



a(fc,r) 



l + 8a(fc,y) In 2 



Pi{k,Yf 



(24) 



with a(k, Y) defined in Eq. ([1^ ). In arriving at Eq. ( p^ ) we had to also include the non-singular terms in a 
in the denominator of Eq. (|23|). This was done for the reasons which will become apparent below, when we 
will explore the convergence of the series we will obtain. 

We can continue construction of our perturbative solution by finding the third order correction N^ik, Y). 
Writing N{k, Y) = Ni{k, Y) + N2{k, Y) + Nsik, Y) we obtain the following equation for Nsik, Y) from Eq. 



dN^ik.Y) 2aNc 



dY 



d 



dink 



aNr 



2Ni{k,Y)N2{k,Y). 



(25) 



The solution to Eq. (25) is given by 



N3(k,Y)^ — 



dy exp 



2aNr 



iy-y)x{- 



d 

dh\k 



[2Ni{k,y)N2{k,y)], 



(26) 
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similarly to Eq. (pT]). Substituting Ni{k,y) from Eq. (15) and N2{k,y) from Eq. (^), performing the 
integration over y and neglecting the two pomeron contribution which arises in the same way as one pomeron 
contribution appeared in Eq. (|2^) we reduce Eq. ( p6|) to the following form 



N3ik,Y) = 2 



dXi d\2 dXs 
(27ri)3 



C\, Ca, Ca, I 



A1+A2+A3 



2[x(-A2)+x(-A3)-x(-A2-A3)] 



1 / nN 

. , ^^ . , ^^ . , ^ , TTT CXp 2 ^ F [x(-Ai) + x(-A2) + xl-As)] ) • (27) 

2[x(-Ai) + x(-A2) +x(-A3) - X(-Ai - A2 - AsjJ V 



After saddle point integrations Eq. (]27|) yields 



iV3(fc,r) = 



a(fc,y) 



2(2 ln2 + !)(! + 8 a(fc,r) ln2; 



.Pi(fc,y)^ 



(28) 



Now we are ready to write down the general expression for the perturbation series of (^^. An obvious 
generalization of Eqs. ( pT| ) and (26) allows us to write 



7v„(fc,y) 



dy exp 



2aNc 



(Y-y)x 



d 



d\nk 



J2 Nm{k,y)Nn-m{k,y) 



(29) 



, m — 1 



As one could see from above the term of order n in the series is proportional to P". Each pomeron splitting 



through the integration over y introduces a denominator of the type present in Eqs. (^3|) and 
series becomes 



The 



iV(fc,r) = ^ /„Pi(fc,r)", 



(30) 



where the coefficients of the series should be determined from the relationship 



71- 

, ^ ^ fm fn — m i 
m— 1 



(31) 



with 



an = < 



^na{k,Y){l + naik,Y) [n2 In 2 - V(l) + V'(n/2)]}"\ 



even n, 



(32) 



i{n2 1n2-V(l)+V'(w/2)}" 



odd 



Eqs. (|30|), (pT] ) and ( p^ provide us with the exact solution of Eq. ( p^ ) in the kinematic region where the 
series of Eq. (^o[) is convergent. Since each term Nn{k, Y) is given by Eq. (|2^) the infinite sum of these terms 
( pO| ) satisfies Eq. (|l^) giving the exact solution to it in the region where it is finite. To obtain Eq. (^) we 
had to assume that na{k, y) <C 1 (as a(fc, Y) <C 1) and expand the x functions in the denominators to the 
first non-singular terms in a(fc, Y). Eq. ( pO| ) represents an expansion of the amplitude in terms of multiple 
hard pomeron exchanges, the nth term in the series corresponding to the n-pomeron exchange contribution. 

Eqs. (^) and (^2|) yield us with the prescription of calculation of an arbitrary term in the series for 
N{k, Y). One could explore the effects of the higher order corrections to the BFKL pomeron arising from 
Eq. (jlO|) by plotting the partial sums of the series of Eq. (|3C|). Defining S'„ — J2m=i /™ ^i"' P^o* S'„ 
as a function of Pi in Fig. ^ for n = 1, . . . , 12. We fixed a to be a constant equal to 0.3 for simplicity. 
The single BFKL pomeron contribution corresponds to the first term in the series, 5*1 — Pi, and is given by 
the straight line in Fig. H, because we plotted Pi as an argument on the horizontal axis. Each partial sum 
Sn gives us some approximation to the exact result for N, and then starting from some value of Pi starts 
diverging. Inclusion of higher orders gives better approximations. That way the thick line formed by the 
overlap of the curves corresponding to different partial sums represents the true value of N, the solution of 
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Eq. (|lO| ) to which the series is converging. One can see that the exact solution shown by that thick line 
lies below the single BFKL pomeron contribution (straight line), deviating from it more and more as the 
energy (and, therefore. Pi) increases. Thus the multiple pomeron contributions start unitarizing the BFKL 
pomeron. An analogous conclusion was reached by Salam in through numerical study of the multiple 
pomeron exchanges in the onium-onium scattering in the framework of the dipole model |l3| . 




FIG. 3. Plot of partial sums (5„) of the series of Eq. (^) as functions of Pi, keeping a = 0.3 for n = 1, . . . , 12. 
The thick line formed by the partial sums corresponds to the true value of N{k, Y). The straight line represents the 
single BFKL pomeron exchange contribution. 

To find the radius of convergence and the high energy asymptotic behavior of the series in Eq. (|30| ) let 
us study the terms corresponding to large values of n. These terms describe the multiple pomeron exchange 
contributions. We would like to explore very large n, such that n^a{k,Y) ^ 1. Then the coefficients a„ 
determining the recursion relations in (^) for the coefficients of the series become 

^ (33) 



n 4 In 2 

Thus, redefining 

ji-i 

.4 ln2, 

we obtain the following recursion relations instead of Eq. ( ^l|) 

n-l 



/" = ( Tli ) /- (34) 



fn ^ - fm fn-m, large 71. (35) 



n 

m—1 



For very large ?i Eq. (35) could be satisfied up to order 1/n corrections with the ansatz 



where r could be some arbitrary constant. Here we have to point out that since Eq. ( |35| ) is written for 
very large n we can not use /i = 1 as a "normalization" condition anymore and that way we have lost the 
information which would have allowed us to fix the value of r. In principle r should depend on a(fc, Y), but 
this function is a much slower varying function of k and Y than Pi(fc, Y), so we assume r to be constant in 
our analysis. If we fix /i = 1 then numerical estimates show that r « 1.4, which allows one to hope that for 
asymptotics of Eq. (^ r is not very large and does not introduce a large correction to our analysis. 

Using the ansatz of Eq. (^) together with Eq. (jsj) in Eq. ( |30| ) we conclude that the constructed 
perturbation series is convergent as long as 
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Piik,Y) 
r4 ln2 



< 1. 



(37) 



Outside the region specified by the condition of Eq. (37) the series in Eq. (pG) is not convergent anymore 



Multiple pomeron exchange contributions corresponding to large values of n become much larger and, con- 



sequently, more important than the one- 
scale Qs by the following condition 



or two-pomeron exchanges. Defining the saturation momentum 



Pi{Qs,Y) 
r4 ln2 



we obtain (for not very large Qg) 



Qs{Y) = A 



exp[(ap ^ 1)Y] 

• 4 In 2 ^UaNcC{3)Y ' 



(38) 



(39) 



Using this quantity we can say that for the transverse momenta of the quark-antiquark pair greater than 
the saturation momentum, k > Qs, the series converges. For the transverse momenta below the saturation 
scale, k < Qs, that is, inside the saturation region, the series diverges ceasing to be a reliable approximation 
to the solution of Eq. ([l^). 

One might wonder whether the condition ln^(fc/A) > aNcY, which we needed to perform the saddle point 
integrations above, would prevent us from exploring the high energy asymptotics with the series of Eq. ( ^ ) 
by putting an upper limit on the possible values of Y for fixed k. Really, this condition provides us with the 
limit on the range of rapidities for the perturbative solution we found 



1 



In" 



A' 



(40) 



However we also have to note that the series (|3^) is convergent as long as fc > QsiY) 
outside of the saturation region. This yields us with another constraint on Y 



-1)Y 



I.e. 



Y < Ysat 



OLp — 1 



In- 



(41) 



Now one can see that the two conditions ( p^ ) and ( [4l| ) are roughly of the same order of magnitude. Thus 
the condition ln^(fc/A) > aN^Y does not impose any additional constraint on the kinematic region of the 
applicability of the series in Eq. (^0|) since the series diverges in the region where the condition is violated 
anyway. 

Now we want to find some approximation to the solution of Eq. ( pO| ) inside the saturation region. The 
saturation region is specified by the condition that k < Qs{Y) and could be achieved either by decreasing 
the transverse momentum or by increasing the center of mass energy of the system. 

First we note that one may find a correction to Eq. (p6[). Rewriting it as 



f - m + e 

Jn — Jn I '^m 

where — (1/f)" for large n and ^ a"/^ ^ 1 for small n, i.e., it satisfies Eq. (|31 
asymptotics given by leading in n behavior of Eq. (|36|). Substituting Eq. (^ ) into Eq. 

to the lowest order in e„ and neglecting for small values of n one obtains 

1 /I " 
e„ ~ - - 
n \r 



(42) 

with the large n 
expanding it 



(43) 



Using Eqs. ([43|), (36) and (|^ in ( pO| ) we could sum up leading in n terms of the series ( |30D obtaining some 
expressions which could be analytically continued outside the region of convergence of the series of Eq. (|30|). 
This way we might obtain an expression for the solution of Eq. (10) inside the saturation region, that is 
for large Pi . While the leading in n part of the series given by Eq. ( |36| ) gives a function which goes to a 
constant for large Pi , the subleading terms given by Eq. (El) yield us with 
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E 



Pi 



r4 ln2 



In 1 



Pi 



r4 ln2 



(44) 



Analytically continuing ln(l + 7^7^72) region of large Pi we obtain for asymptotically large energies 

course we do not know whether the higher order corrections in 1/n would not give us some 
other function of Pi which would be large for large Pi possibly interfering with our In Pi ansatz. Therefore 
In Pi should be considered only as a a possible high energy asymptotics of N and the above argument should 
not be taken as a proof of this statement. 

To justify that the above ansatz of iV = In ^. ^^^^ ^ = In ^ is the solution of Eq. (^0|) at asymptotically 
high energies we first note that 



N{k,Y) = In 



M 



(45) 



with M some large momentum scale (k < M) satisfies the equation 



2X - 



d 



ainfc 



N{k,Y) -N^{k,Y) = 0. 



(46) 



Really, rewriting 



M 



dX 1 
A2 



(47) 



where the contour C is a circle of radius less than one centered around A = 0, one can see that since x(— A) 
also has a pole at A = then 



2x 



d\nk J k k 



(48) 



We can assume that in the saturation region N{k, Y) ceases to depend on Y. Then Eq. (10) reduces to Eq. 
(^6|). To find the most general solution of Eq. (46) we represent TV as a series 



iv(fc,r) 



n.ni—O 



M 



(49) 



and substitute it into Eq. (^6|) in order to find the unknown coefficients Cnm- Using the technique outlined 



above for In ^ we can show that 



2x - 



dink 



— V In" — 

M I ^ k 



2 /'AV* lr,™+l M. 
m+1 \ m) k ' 

2x(-n)(A)"ln"f, 



even n. 



odd 



(50) 



With the help of Eq. (50) one can show that in order for the series of Eq. ( [49|) to satisfy Eq. ( [46[ ) all of 
its coefficients must be zero with the exception of cqi = 1. Thus we proved that the most general energy 
independent solution of Eq. ( p^ is given by Eq. (^5|). 

At asymptotically high energies there is only one large momentum scale in the scattering of the quark- 
antiquark pair on a hadron or nucleus — the saturation scale Qs ■ This scale is also s ugg ested by the ansatz 
of Eq. (^). Therefore we should put M = Qg and use Eq. ( ^ ) as a solution of Eq. ( |l(]| ) in the region when 
k < Qs- The derivative of this new solution with respect to rapidity is of the order of dln(Qs/k)/dY ^ aNc 
and could be neglected compared to any of the terms on the right hand side of Eq. (^^, each of them 
being of the order of aNc{aNcY)'^ . This could be done since aNcY ^ 1. Therefore we conclude that in the 
saturation region 



N{k,Y) 



, Qs 



Y > 



1 



ap — 1 



ln4 



(51) 
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Now that we know the solution for N{k, Y) outside of the saturation region given by Eqs. (^o|), and 
( ^ ) and the high energy asymptotics inside the saturation region given by Eq. (51), one might want to 
construct the solution in coordinate space {N{x±,Y)) using the Fourier transformation of Eq. (pa|). To do 
that one would have to integrate over all values of the transverse momentum k from to oo. If x± < 1/Qs 
then the transverse momentum, which in Eq. ( ^ ) is effectively cut off by the Bessel function and , therefore, 
varies approximately from to l/x± > Qs, will go through a range of values both above and below Qs- For 
the first case one would have to use the perturbation series which we constructed as N, while for the second 
case one might use the asymptotic value of N given by Eq. (|5^) , although a perturbative expansion around 
it would be more accurate. The result would be a complicated combination of special function and we are 
not going to list it here. 

If x± > 1/Qs then the transverse momentum in Eq. ( pa| ) simply gets cut off by l/x± < Qs and always 
stays within the saturation region. Thus one could just perform a Fourier transformation of formula ( pT| ) 
using Eq. (|9a|). The result is 

N{x^,Y) « 1, Y> In \. (52) 

This corresponds to the blackness of the total cross section of the quark-antiquark pair on the target nucleus, 
since it is given by 

alf^ ^ 2 J (fhoN{x^M. Y) « 27ri?2 (53) 

for N = \ and for a cylindrical nucleus of radius R as described above. That way the total cross section is 
independent of energy at asymptotically high energies and is completely unitary. It reaches its "geometrical" 
limit (BSh which could be predicted even from quantum mechanics. 



1 




FIG. 4. Plot of the qualitative behavior of the forward qq pair scattering amplitude on a nucleus, N{x_i_,Y) , as a 
function of x±Qs{Y) (one BFKL pomeron exchange contribution). At small values of x±Qs the structure function is 
linear in x±Qs, reproducing BFKL pomeron. As energy increases the forward amplitude A'' saturates to a constant. 

That way we have shown that N{x±,Y) given by the solution of Eq. (|l^) behaves like a single BFKL 
pomeron exchange contribution at moderately high energies (Y ^ 1 /a) , which follows from the perturbation 
series we have constructed, and as energy increases to very high quantities (Y > In ^), saturates to 

a constant independent of energy. A qualitative sketch of N{x±,Y) as a function of x±Qs(Y) (one pomeron 
contribution) is shown in Fig. 



III. ASYMPTOTIC BEHAVIOR OF THE F2 STRUCTURE FUNCTION 



For a large cylindrical nucleus the F2 structure function is given by 
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(fxoi dz $(xoi, z) N{x(^i,Y), 



which follows from Eq. (0). We can rewrite Eq. (Is^) in momentum space using Eq. (pa|) 



where 



$(fc,z) 



2\ Q^R^ 



EM 



(fkdz^{k,z)N{k,Y), 



(2^)2 



Xqi ^{xoi,z). 



Employing Eqs. ( |5aD and (5b) we obtain 



(54) 



(55) 



(56) 



^k, z) = ^ { [z^ + (1 - zn ^ U 3, 2.5, I + 



+4Q^z^(l - zY 



15afe 



5 a' 2F1 I 2, 2, 2.5, --^ ]-2k\Fi{ 3, 3, 3.5, — ^ 



4a2 



4a2 



(57) 



where 2-^1 is a hypergeometric function. 

Eq. (^5]) shows that in order to write down an expression for the F2 structure function one has to know 
N{k, Y) in the areas above and below saturation. The situation is similar to the Fourier transformation 
of the previous section. The wave function $(fc, z) given by Eq. (pTtj becomes very small for transverse 



momenta k > a = Qy^z{l — z), effectively providing an upper cutoff on the k integration. Thus ii Q > Qs 
the integral over k includes k < Qs and k > Qs and to calculate it we have to use the solution for N both 
inside and outside of the saturation region. We can again propose the approximation outlined in Sect. II 
for calculation of F2, which consists of using perturbation series of Eq. ( |30| ) outside saturation region and 
approximating the solution inside that region by Eq. (^). The result would also include a complicated list 
of special functions which we are not going to list here. We have to admit that numerical solution of Eq. 
(^ could give a more precise results for F2. Here we can just point out that at moderately high energies the 
behavior of the F2 structure function will be dominated by single BFKL pomeron exchange, providing us 
with the following formula which could be obtained by substituting the first term in the series of Eq. ( ^ ) 
into Eq. ^) 



F2ix,Q^)^^S^N,QA 
255 



C_i exp[(ap-l)r] 



4 1n2 y/UaNX{3)Y ' 



(58) 



where we have neglected the transverse momentum diffusion term in the exponent of Eq. (|l8|) and S± — ttR^ . 
Eq. ( |5^ is valid only for moderately high energies, when the saturation momentum Qs is small and we can 
neglect the k integration in Eq. ( |55| ) below that scale. 

When Q < Qs or the energy gets very high to sufficiently increase Qs the integral over k in Eq. ( ^5| ) is 
limited to fc < Q < Qs, so that we can use Eq. (pl|) for N{k, Y). The result yields 



F2{x,Q^) 



Q^R'N^ 
37r2 



In 



Qs 
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Q <Qs 



(59) 



As energy becomes very large, corresponding to Y ^ ap-i t?' ^ becomes approximately equal to 
(ap - 1)Y. Therefore 



F2{x,Q^ 



Q^R'Nc 
37r2 



{ap - i)r. 



Y > 



ap — 1 



■In- 



(60) 



We conclude that F2 ~ Ins at asymptotically high energies. This conclusion may seem to be a little unusual, 
so we will reproduce the same result using more conventional coordinate space language. As was shown in 
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the previous section the forward amphtude of the qq pair scattering on a nucleus is N{x±,Y) = 1, when 
x± > 1/Qs- When x± < l/Qs one could argue that the amplitude is roughly proportional to some positive 
power of x±, at least N{x±,Y) ^ x±, as given by one pomeron exchange. Thus, for Q < Qs we can neglect 
the part of the integral with xqi < 1/Qs in Eq. (p4|), since the wave function in Eq. (jsj) behaves like 
and is not singular enough to make that portion of the integral significant. That way, expanding the 

modified Bessel functions in $ of Eq. (|5a|) and integrating over z we obtain 



F2(a;,Q)« — ^ — / . (61) 



37r2 



The upper cutoff on the xqi integration in Eq. ( |6l| ) is provided by the fact that the modified Bessel functions 
fall off exponentially at the large values of the argument. Integrating over xqi in Eq. ( |6l| ) we arrive at Eq. 
®. 

That way we have shown that the F2 structure function given by the solution of Eq. (g9|) is proportional 
to the single BFKL pomeron exchange contribution ( ^8|) at moderately high energies with F ~ 1/a, and, as 
energy increases the F2 structure function unitarizes, becoming linearly proportional to Ins (^0|). We stress 
again that even though the cross section of the quark-antiquark pair of a fixed transverse size x± scattering 
on a nucleus saturates to a constant at large Y, the integral over xj_ makes F2 depend on Y, such that 
F2 ^Y dX very large Y . 



IV. APPLICATION TO DIFFRACTION 



In Jp5[ a diffractive deep inelastic scattering process was considered, an amplitude of which is shown in 
Fig. pTA virtual photon scatters on a nucleus elastically, leaving it intact. Two quarks in the final state in 
Fig. H may become a pair of jets or form a vector meson. These type of processes were called quasi-elastic 
virtual photoproduction in p^ ]. 




••• C3 

FIG. 5. Quasi-elastic deep inelastic scattering as pictured in [25]. 

The interactions between the quark-antiquark pair and the nucleus were taken in the quasi-classical 
approximation, similar to p^ . Each nucleon in [ p5[ would interact with the qq pair by a two gluon exchange. 
Strictly speaking this approximation corresponds to not very high energies, when the logarithms of s resulting 
from QCD evolution are not large enough to become important. That is why the approximation is called 
quasi-classical. In [ p5[ it provided us with a simple Glauber typ e expression for the diffractive, or, more 
correctly, quasi-elastic structure function F2 [cf. Eq. (20) in |2^ ]] 

J^2''(^,Q') = ^^^ / "^'^"i j\z{a^Kl{xo^a)[z^ + {\-zf]+AQ^z^{\-zfKl{xo^a)] 



X72(xoi,bo), (62) 

where we again assumed that the quarks are massless and have only one flavor. 7(xoi,bo) is the Glauber 
propagator of the quark-antiquark pair, which for the case of a cylindrical nucleus is given by Eq. (0). One 



15 



can obtain Eq. ( |62|) from Eq. (y) by first noting that the quasi-classical limit corresponds to putting Y = 
(no evolution in rapidity) in it and iV(xoi,bo,y = 0) = 7(xoi,bo). Recalling that, as was argued in p5[ |, 
one can rewrite the total and elastic cross sections for DIS on a nucleus in terms of the S'-matrix at a given 
impact parameter of the collision, S{b), yielding 



atot ^2 d?b[l- S{h)] 



(63a) 



n,i = / (fb [1 - S{b)f 



(63b) 



and associating F2 with atot and with aei , we can easily derive Eq. (|6^) from Eq. (^) . 

Now let us imagine our quasi-elastic process at very high energies. Then, strictly speaking, one has to 
include the evolution in rapidity in the structure function F2 ■ We still do not want to have particles being 
produced in the central rapidity region. One could model the interaction by a single pomeron exchange be- 
tween the quark-antiquark pair and the nucleus, similarly to what was done previously in and references 
mentioned there for elastic scattering. Here we are going to propose inclusion of multiple pomeron exchanges 
in . Since the above formalism allows us to sum up pomeron fan diagrams in the leading logarithmic 
approximation, we can just insert this evolution in the amplitude of the quasi-elastic process. The diagram 
of the process is depicted in Fig. |6| The diffractive rapidity gap between the nucleus in the final state and 
the remains of the qq pair results not only from a single pomeron exchange, but also from multiple pomeron 
exchanges. Still nothing is produced in the central rapidity region. 




FIG. 6. Quasi-elastic scattering including multiple pomeron exchanges, as envisioned in our approach. 

The corresponding expression for the diffractive structure function F^ is easy to derive. We know the 
expression for the total F2 structure function in Eq. (^, which includes all the discussed multiple pomeron 
exchanges. One can easily argue along the lines outlined in deriving Eq. (|62|), that for the case of multiple 
pomeron exchanges we a gain have a relationship between the total and elastic cross sections similar to the 
one given by Eqs. (|3|) an d (^ ). That way we conclude that to obtain a formula for F^ one has to 
substitute 7(xoi,bo) in Eq. (p2) with A^(xoi, bo, F). The result is 



7r(27r)3 



d'^bo d^xoi / dz {a^Kf{xoia)[z^ + (1 - zf] + AQ'^z^{l - zfK'^{xoia)} 



(64) 
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where 7V(xoi,bo,y) is, in the large Nc hmit, a solution of Eq. (||). Eq. (^) is our master formula for the 
diffractive (quasi-elastic) structure function. When Y = (no evolution) it reduces to Eq. (|6^). 

Using the results of Sect. II we can derive the asymptotics of the diffractive structure function. At 
moderately high energies the one pomeron exchange contribution dominates the diffractive amplitude of 
Fig. 0. Again, similar to the previous section we assume that since at y 1/a the saturation scale Qs is 



small one can neglect the k integration below Qg. Performing a Fourier transformation (9a) of the first term 



in the series of Eq. ( pOD and substituting the result into Eq. (|64| ) we obtain for a cylindrical nucleus 

^"(-'Q) = ^ (4 1n> lWc(3)r ^' ^^a' ^''^ 

where we had to cut off the z integral in Eq. ( ]6^ ) by e^^ and we again neglected the transverse momentum 
diffusion term in Eq. (|l^ . The diffractive structure function of Eq. (|65| ) behaves like a square of a single 
pomeron exchange contribution, similar to what was found for elastic parton scattering in p9|. Even though 



P'p ~ Ql(X) ™ Eq. (|65| ) is growing with energy faster than F2 ^ QsiY) in Eq. (58), unitarity is not 
violated, since both expressions are valid in the region of relatively small Qs corresponding to Y 1/a. 
Thus at moderate energies <F2, as < A. 

We also note that formula ( |65| ) was derived assuming that the photon's virtuality is not very large. 
Therefore, the fact that in Eq. (|65| ) has no powers of Q in the numerator, whereas the corresponding 
F2 given by Eq. ( p8| ) is proportional to Q, does not imply that FP is not a leading twist effect, violating 
the theoretical predictions ||2^,0 and experimental observations |^ . For the case of large Q one should 
perform a different saddle point approximation for the terms in the perturbative series given by Eqs. (p^), 
(p^), etc., [see Q], getting a different series solving Eq. ( px| ) outside of the saturation region. This should 
result in concluding that F^ given by the first term in that series is not suppressed by any powers of Q 
compared to -Pj- However, doing that is beyond the scope of this paper. 

At very high energies, when Qs{y) is large enough so that Q < Qs{Y), the quasi-elastic structure function 
F^ saturates. Similarly to Sect. Ill we can derive the high energy asymptotic behavior of Eq. (|64|). Putting 
N{xoi,Y) = 1 for xqi > 1/Qs and imposing the cutoffs on the xoi integration in Eq. (^J) analogously to 
what was done in deriving Eq. (|6l|) we deduce 



i^2^(x,Q^)«^^^(ap-l)y, Y>^—ln^. (66) 



That way the F^ also becomes linearly proportional to the logarithm of center of mass energy energy at 
very high energies, similarly to the total structure function i^2- Comparing Eq. (^6|) with the high energy 
behavior of the F2 structure function given by Eq. ( |60| ) we conclude that F^ is still smaller than F2 and 
unitarity is again not violated. One can see that at this asymptotically high energy F2 = 2F^, which 
corresponds to the prediction for the black cross sections resulting from quantum mechanics. 



V. CONCLUSIONS 

In this paper we have found a solution of the evolution equation (|^) which resums all pomeron fan diagrams 
in the leading logarithmic limit for large that was proposed in [Q. The solution consists of two parts. 
Outside of the saturation region we have constructed a perturbative series in momentum space given by 
Eqs. (^), ( |3l| ) and (^2|). We have demonstrated that the series tends to unitarize the single BFKL pomeron 
exchange contribution (see Fig. |3|) . In deriving formula which led to Eqs. (|l|) and (||), we had to 

assume that DIS cross sections are slowly varying functions of the impact parameter, which corresponds to 
scattering on a large nucleus. Only in the large nucleus case we could argue that the fan diagrams dominate 
the total DIS cross sections. In order to perform the saddle point approximation in Eqs. (p^), ([22[), ( p7| ) and 
for the series ( |30| ) to be convergent we had to assume that ln(Q/A) <C aNcY and Q > (YTT^restricting 
the possible values of Q to a certain not very broad range. Nevertheless the abovementioned assumptions 
did not change the qualitative behavior of the obtained result and the conclusions could be easily extended 
to large Q by performing a different saddle point approximation. 

Inside the saturation region, for Q < Qs, we found an approximate solution given by Eq. (^) correspond- 
ing to the saturation of the total scattering cross section of a quark-antiquark pair of a given transverse 
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size on a nucleus to a constant independent of energy That way we have shown that the coordinate 

space forward ampHtude of the qq scatterina on a nucleus or hadron grows like a BFKL pomcron exchange 
contribution at moderately high energies (18), and, as energy gets very large it saturates to a constant (|5^). 
A qualitative plot of N{x^,Y) is presented in Fig. ^. 

We have to admit that in order to construct a solution of Eq. (^) in coordinate space, N{x±, Y), and the 
corresponding structure function F2 one has to have a better knowledge of the momentum space solution 
inside the saturation region. A numerical solution of Eq. would probably be very helpful in determining 
exact values of N{x±,Y) and F2 at intermediately large rapidities. 

We have shown that the F2 structure function given by the solution of Eq. grows as a power of energy 
at moderately high energies, corresponding to rapidities of the order of Y ^ 1/a. This behavior, presented 
in Eq. (|58|), reproduces the usual single BFKL pomeron exchange contribution to the structure function. 
As energy grows very high, corresponding to the rapidities of the order of F> [l/{ap — 1)] ln(l/Q!^), the F2 
structure functions unitarizes, becoming linearly proportional to Ins, which satisfies Froissart bound. This is 
shown in Eq. @. Thus we have shown that the BFKL pomcron on the nucleus is unitarized by summation 
of pomeron fan diagrams. 

We also note that we derived an expression for the quasi-elastic (diffractive) structure function F^ pre- 
sented in Eq. (^), which also includes all multiple pomeron exchanges. We have shown that at moderately 
high energies the obtained expression behaves like a square of the BFKL pomeron contribution [Eq. (|65|)] 
and at very high energies it saturates and becomes a linear function of Ins [Eq. (|66|)], being equal to a half 
of the structure function F2, as expected for black cross sections. 

In Q it was argued that the dipole model provides us with the techniques that allow one to resum pomeron 
fan diagrams at any subleading logarithmic order in the large Nc limit. That is, if one calculates the next-to- 
leading order kernel in the dipole model, the resulting equation for the generating functional Z [see [pl]-p^] 
would provide us with the equation resumming all NLO pomcron fan diagrams. The resulting equation 
would be, probably, more complicated than our Eq. (H). It would have a much more sophisticated kernel 
[0,|j and would have cubic terms in N on the right hand side. However, the analysis presented in this paper 
allows one to hope that this resummation will cure some of the problems of the recently calculated NLO 
BFKL kernel [^ . Resummation of one loop running coupling corrections, which are a part of the NLO BFKL 
kernel, brings a factor of e*^" ^ in the one pomeron exchange amplitude, with the constant c ~ 5 for three 
flavors |9j . Since it introduces an even faster growth with energy than the leading order BFKL pomeron, this 
factor seems to complicate the problem of unitarization of the pomeron. However, the inclusion if this factor 
of e'^" ^ in our definition of Pi in Eq. (|l8|) would not change the perturbation series of Eq. ( pO| ) and the 
quark-antiquark scattering amplitude would still be unitarized by multiple pomeron exchanges. Of course 
a careful inclusion of higher order corrections would lead to much more complicated results. Nevertheless 
we may hope that inclusion of one loop running coupling corrections would just modify the one pomeron 
exchange contribution and the multiple pomeron exchanges would still unitarize the total DIS cross section. 
However, a rigorous proof of this statement is beyond our goals here. 

Finally we point out that one can use the developed technique to construct a gluon distribution function 
xG{x, Q^) of the nucleus including multiple hard pomeron exchanges similarly to the way we have constructed 
the F2 structure function. To do that one has to consider scattering of the "current" j = — (l/4)F^^,F^^ on 
the nucleus, similarly to |l^,^,^. Then, in the large Nc limit, it should be possible to derive an equation 
for the generating functional of the dipole wave function of the current j, which, at the end, would provide 
us with an expression for the gluon distribution function xG{x,Q^). 

After this paper was finished the author has learned about another effort to solve Eq. (||) which was 
carried out in |33l. 
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APPENDIX A: 



Here we are going to show how the Fourier transformation of Eqs. ( Pa|) and 
(0). First we note that 



reduces Eq. m) to Eq. 



d X2 = "^n xi2 xo2 dxi2 dxo2 / dkk Jo{kxoi) Joikxo2) Joikxi2), 



(Al) 



which would allo w u s to integrate independently over Xq2 and xi2 in the kernel of Eq. (g). Using the 
Jacobian of Eq. (Al) in the second (quadratic) term on the right hand side of Eq. (^), multiplying both 
sides of Eq. m) by Jo(fc2:oi)/a;oi, integrating over a;oi from to oo and employing 



dx X Jo(kx)Jo(k'x) — —S(k~k') 
k 



(A2) 



we obtain 



dNjk.Y) aN, dxoi ^ , 



dY 



a;oi 



d X2 



^02^12 



2-11:5 (xoi - X02) In 



A^(xo2,r) 



-N^{k,Y). 



(A3) 



The first (linear in N) term on the right hand side of Eq. (A3) can be rewritten employing Eqs. (Al), ( |9a| ) 
and (9b) as 



2aNc 



dk'k'N{k',Y) / dxoi xqi Jo{kxoi) Jo{k' xqi) 



dxi2 , , , xoi 
Jo(fc X12 ) - In 

Xl2 p 



(A4) 



Using the techniques outlined in Eqs. (32) and (33) of ||ll| the integral in the parenthesis of Eq. (A4) can 
be done, yielding 



Jo(fc 2^12) = V'(l) - In 

Xl2 2 



(A5) 



Substituting this back into Eq. (A4) we end up with 
2aN, 



dk'k'N{k',Y) / dxoiXoiJQ{kxoi)Jo{k'xoi) ^(1) - In 



k'xQi 



Expanding Jo(fca:oi) in a Taylor series we rewrite Eq. (A6) as 



2aNc 



dk'k' N{k',Y) ^ 



(-1)" 
(m!)2 



^(1) 



d 



d{2m) 



dxoi xoiJo{k'xoi)xl'^ 



Performing the integration over xqi in Eq. (A?) we obtain [ see Eq. (33) in pTj] 



2aNr 



°^ dk' ~ , , ^ (-1) 
-N{k',Y)Y,^ ' 



k' 



— ' (rn!)2 

m=Q ^ ' 



d 



d{2m) 



2m 



r(m + 1) 
r(-m) 



After differentiating over m and employing the fact that 



E 



(-1)'^ 



to! 



(A6) 



(A7) 



(A8) 



(A9) 



which is a direct consequence of Eq. (A2), we can reduce Eq. (A8) to 
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N{k\Y)J2'—^2x{-2m) 



r(-m)' 



(AlO) 



with X defined by Eq. ( [Tl| ) above. Rewriting x(— 2m) in Eq. (AlO) as x{~9/dhik) and employing Eq. (A9) 
to perform the summation, which makes the integration over k' trivial, we obtain the following expression 
for the first term on the right hand side of Eq. (A3) 



2aNr 



d 



d\nk 



N{k,Y). 



(All) 



That way Eq.(^) reduces to Eq. M), as desired. 
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